Entanglement induced in spin-i particles by a spin-chain near its critical points 
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A relation between entanglement and criticality of spin chains is established. The entanglement 
we exploit is shared between auxiliary particles, which are isolated from each other, but are coupled 
to the same critical spin- 1/2 chain. We analytically evaluate the reduced density matrix, and 
numerically show the entanglement of the auxiliary particles in the proximity of the critical points 
of the spin chain. We find that the entanglement induced by the spin-chain may reach one, and 
it can signal very well the critical points of the chain. A physical understanding and experimental 
realization with trapped ions are presented. 
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Quantum entanglement lies at the heart of the difFer- 
ence between the quantum and classical multi-particle 
world, and can be treated as a useful resource in various 
tasks such as cryptography, quantum computation and 
teleportation [^. Quantum phase transitions [2] are tran- 
sitions between quantitatively distinct phases of quantum 
many-body systems, driven solely by quantum fluctua- 
tions. In the past decade, a great effort has been devoted 
to understand the relations between entanglement and 
quantum phase transitions [3, 0, S S, 0, HQ- In fact, 
it is natural to associate the quantum phase transition 
and entanglement once correlations are behind both of 
them. By sharing this point of view, one anticipates that 
entanglement induced by a quantum critical many-body 
system will furnish a dramatic signature of quantum crit- 
ical points for the many-body system. 

On the other hand, we usually think of environment 
that surrounds quantum system as a source of decoher- 
ence. Recently researchers have started to investigate the 
positive effects [13, [nl, H [S Ei, El, H El of environ- 
ment, for example, environment assisted information pro- 
cessing and environment induced entanglement. These 
investigations pave a new way to engineer mechanisms 
of preventing, minimizing or using the impact of envi- 
ronment in quantum information processing. In those 
works, however, the environment was modelled as a set 
of independent quantum systems, i.e., correlations among 
particles in the environment were ignored. An interesting 
open question is whether the correlation among environ- 
mental particles can affect the entanglement induced in 
a bipartite system that couples to it. 

In this paper, we show how to exploit entanglement in 
auxiliary particles induced by a quantum critical many- 
body system as an essential tool to reveal quantum phe- 
nomena in the many-body quantum system. Indeed, 
quantum phase transitions are accompanied by a quali- 
tative change in the nature of classical correlations, such 
drastic changes in the properties of ground states are 
often due to the coUectiveness/randomness of interpar- 
ticle couplings which are possibly reflected in entangle- 
ment between systems that couple to it. Here we adopt a 



spin-chain system described by the one-dimensional spin- 
■i XY model as the many-body system. Another pair of 
spin-i systems that couple to the spin chain would act as 
the auxiliary particles. We observe that the entanglement 
of the auxiliary particles changes sharply in the proximity 
of quantum phase transition. This change can be traced 
down to the presence of coUectiveness in the dominant 
couplings of the auxiliary particles to the chain, then it 
reflects the critical properties of the many-body system. 
This observation offers a new tool to study quantum crit- 
ical phenomena, in particular for more general systems, 
where analytical solutions might not be available. A pos- 
sible realization of this scheme is proposed with trapped 
ions. It utilizes off-resonant standing- waves driven ions 
in traps [ll,E3| to simulate the one-dimensional XY spin- 
chain. The proposal also could be realized with ultracold 
atoms superposed by optical lattice (23|. This makes the 
study appealing for experimental test. 

Consider a bipartite system consisting of two spin-i 
particles a and 6, and a quantum many-body system de- 
scribed by the one dimensional spin-i XY model. The 
system Hamiltonian Hs and the Hamiltonian Hb of the 
chain read, 
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where N is the number of sites, erf (a — x,y,z) are the 
pauli matrices, and 7 is the anisotropy parameter. The 
periodic boundary condition (Tat+i — ai is assumed for 
the spin chain. Suppose the coupling of the auxiliary 
particles to chain take the form 
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where g and h denote coupling constants. Clearly, 
[Hs, Hi] — 0, which implies that the energy of the auxil- 
iary particles is conserved, but the coherence may not be 
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conserved depending on the detail of the system-chain 
coupHng. This leads to the following form of the time 
evolution operator U{t), U{t) = J2i,j=o,iUij{t)\ij){ij\, 
where = |i)a «) \j)b, and \i)a{\ijb,i = 0,1) repre- 
sent the eigenstates of (J^{<J^)- It is easy to show that 
U^j{t){i,j = 0,1) satisfy ihpj,j{t) = H.^U^it) with 

= -E;ii(^^f^f+i + ^^/Vf+i + ^^l). where 
= A (-l)*+i2g + (-l)J+i2/i, i,j = 0,1. If the 
auxiliary particles are initially in state the dynam- 
ics and statistical properties of the spin chain would 
be govern by Hij, it takes the same form as Hb but 
with modified field strengths Ay . The Hamiltonian Hij 
can be diagonalized by a standard procedure [IJl to be 

= 'Ek'^i3Avl,kV^j,k - 5), where mj.k{4j,k) are the 
annihilation (creation) operators of the fermionic modes 

with frequency ujij^k 

A,y), k = - 



(cos 2^ - A,y), k = -N/2, -N/2 + 1, iV/2 - 1. The 
fermionic operator rjij^t^ was defined by the Bogoliubov 
transformation as, rjij^k = rffc cos — idlj,sin^^, 
where dk — -^^^ X]; '^iexp(— i27rZfc/A^), and the mix- 
ing angle was defined by cos 
Fermionic operators a; were connected with the spin op- 
erators by the Jordan- Wigner transformation via ai = 



(n.m<i'^m)i^f + ^f^f )/2- The operators rjij^k parame- 
terized by i and j clearly do not commute with each 
other, this will leads to entanglement in the auxiliary 
particles as shown later on. Before going on to calcu- 
late the reduced density matrix, we present a discussion 
on the diagonalization of Hij. For a chain with peri- 
odic boundary condition, i.e., cti — a^, boundary terms 
Hboun ~ [(ajv^i + 7ajyQi) +_h.c.]\exp(iTrM) + 1] have to 
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23j, because we are interested 



be taken into acount U 
work with Hboun ignored 
in finding a link between the criticality of the chain and 
the entanglement in the auxiliary particles. 

Having given an initial product (separable) state of the 
total system, |V'(0)) = \M0)) «) \M0)) «)|'/'s(0)), we 
can obtain the reduced density matrix for the auxiliary 
particles as pafc(t) = TrB[C/(i)l'/'(0))(^(0)|C/t(i)], it may 
be formally written in the form 



Pab{t) = PiJ-:mn{t)\ij){r 
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A straightforward but somewhat tedious calculation 
shows that 
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where cos 9k = 



cos(27rfc/jV)-A 



To derive 
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this result, the spin chain was assumed to be initially 
in the ground state of Hb- Discussions on Eq.(|4]) are 
in order. For g = h, one has Ay = A when i — and 
j = 1, or i = 1 and j = 0. Hence, |01) and |10) expand 
a decoherence free subspace. Subsequently, the entangle- 
ment of state in this subspace remains unchanged due to 
^ij;mn{t) — 1- The situatiou changes if h ^ g, where the 
decoherence free subspace does not exist. The entangle- 
ment shared between the auxiliary particles evolves with 
time according to Eq. (g]) in this situation. 

With these expressions, we now turn to study entan- 
glement shared between the auxiliary particles a and b 
in state ([3]). To be specific, we choose \(j)a{0)) \4>b{0)) 
= l/\/2(|0)a-h|l)a (gil/V2(|0)6-h|l)f, as the initial State of 
the auxiliary particles, while the spin chain is assumed to 
be in the ground state of Hamiltonian Hb- The entan- 
glement measured by the Wootters concurrence can be 
calculated and the numerical result was shown in figures 



1-4. The regions of criticality appear when the ground 
and first excited states become degenerate. We first fo- 
cus on the criticality in the XX model. The XX model, 
which corresponds to 7 = 0, has a criticality region along 
the line between A = 1 and A = — 1. The criticality is 
refiected in the entanglement of the auxiliary particles, 
which appear in figure 1 and 2. Figure 1 shows the Woot- 
ters concurrence as a function of time and the anisotropy 
parameter 7, it is clear that the Wootters concurrence 
has a sharp change in the limit 7^0. This result can be 
understood by considering the value of Oij^k, which take 
or TT depending on the sign of cos(27rfc/7V) — A^ in this 
limit. In either case |ry;mn(OI — 1; which indicates that 
the norm of any element of the reduced density matrix 
remains unchanged. Figure 2 shows the entanglement 
of the auxiliary particles near the criticality region with 
7 = and A = ±1. The entanglement changes sharply 
along the line of A = ±1 in the time-A plane. This is a 
refiection of the critical phenomena in the entanglement 
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FIG. 3: (color online)Entanglement(measured by Wootters 
concurrence) induced by the transverse Ising spin-i chain, 
which may be obtained from the XY model by setting 7 = 
1. The numerical calculation was performed with A'^ = 1502 
sites, and g = h = 0.05. 



FIG. 1; (color online)Wootters concurrence of auxiliary par- 
ticles as a function of time and the anisotropy parameter 7. 
The figure was plotted for A'^ = 1502 sites, and g = h — 0.05. 
Figure l-(a),(b),(c)and (d) correspond to A = 5, A = 2, A = 1 
and A — 0.8, respectively. 
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FIG. 2: (color online) Wootters concurrence versus time and 
A with fixed 7 — 10^^ (close to zero). This corresponds to 
the XX model, and the critical points with A = ±1 are clearly 
shown. The figure on the left is a contour plot for the right. 
The other parameters chosen are A*' — 1502, and g = h = 0.05. 



of the auxiliary particles. 

Next we turn to study the criticality in the trans- 
verse Ising model. The structure of the transverse model 
ground state changes dramatically as the parameter A 
is varied. The dependence of the entanglement on A is 
quite comphcated[3|- Here we present analyses for the 
A — !■ 00, A = 1 and A = limits. In the A — > 00 limit 
the ground state of the chain approaches a product of 
spins pointing the positive z direction. The mixing an- 
gle 6ij,k and 9k tend to — tt in this limit. This leads to 
^ij;mn{t) — 1, which indicates that the initial state does 
not evolve with time, i.e., the auxiliary particles remain 
in separable states. The A = limit is fundamentally 
different from the A — *■ oo limit because the correspond- 
ing ground state is doubly degenerate under the global 
spin flip by Y[i=i ■ This symmetry breaks at |A| = 1 
and the chain develops a nonzero magnetization (cr^) ^ 
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FIG. 4: (color online) Wootters concurrence at time t — 0.75. 
The locations of peaks on the A axis do not change with time. 
Various parameters are A'^ = 1502, g = h = 0.05. 



which grows as A is decreased. This dramatic change in 
the ground state of the spin chain can be found in the en- 
tanglement of auxiliary particles in figure 3. In fact, the 
ground state of the XY models is v ery complicated with 
many different regimes of behavior 2J, 2^ . With what- 
ever 7, there is a sharp change in the entanglement across 
the line |A| = 1 (figure 4). This signals the change in the 
ground state of the spin chain from paramagnetic phase 
to the others. Although the entanglement in both cases 
with 7 = and 7=1 shares the same properties along 
the line of |A| = 1, i.e., it changes dramatically across 
this critical region, the two case are quite different in its 
nature. The two auxiliary particles never evolve in the 
case of 7 = 0, while the two will eventually evolve into 
the pointer states in another case(except for A —>■ oo), 
hence no entanglement share between the auxiliary par- 
ticles in the later case in the t — > oo and N oo limit. 



The connection between the entanglement in auxil- 
iary particles and the criticality of spin chain may be 
understood as singularity in one or more pij-mn(t) (the 
elements of the reduced density matrix) at the critical 
points. To show this point, we recall that 



Pij;mn{t) = - ((/)_B (0) | e" 



-iHijt^iHkit 



'*I</'b(0)), (5) 
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with the same notations and initial states given before 
Eq.([3])- Notice that |(/)b(0)) was taken as the ground 
state \0)b of Hb, but it is not an eigenstates of Hij with 



Aij 7^ A. Expanding \0)b in terms of eigenstate 
Hij, (a = 0, 1, ) 



\0)i 



of 



(6) 



one can easily prove that |c^_^q| ^ h,g when g <C 1 and 
/i <C 1. This is exactly the case under our consideration. 
Where |0)^ denotes the ground state of Hij. Therefore, 



Pij;mn{t) ^ ^ Cq-* [Cq 



l{0\0)] 



(7) 



up to first order in g and h. Here fiy stands for the 
ground state energy of Hij . As observed in [2^ , a sudden 
drop in g (0|0)^ can signal the regions of criticality in the 
spin chain. Therefore the entanglement as a function of 
Pij;mn (t) can signal the criticality of the chain. 

This discussion, apart from its theoretical interests, of- 
fers a possible experimental method to study critical phe- 
nomena without the need to identify the state of the sys- 
tem, in particular in the presence of degeneracy. The XY 
model can be realized with tra ppe d ions under the action 
of off-resonant standing waves [l9j as follows. Consider a 
system of trapped ions, whose physical implementation 
corresponds to Coulomb chains in Penning traps or an ar- 
ray of ion microtraps. The Coulomb repulsion together 
with the trapped potential and ions' motion yield a set of 
collective vibrational modes, these collective modes can 
be coupled to the internal states of ions by conditional 
forces. The directions of the conditional forces [l9l| would 
determine the effective couplings between the effective 
spins. To simulate the XY model, we apply two condi- 
tional forces in both directions x and y. The effective 
coupling is proportional to where r represents the 

relative distance between ions. Although the couplings 
are distance dependent, the critical properties are shown 
to be similar to these ideal models 1^ 27 1 considered in 
this paper. The coupling of the auxiliary particles, which 
share another trapping with the ions in the chain, to the 
chain can be simulated in the same way, where a condi- 
tional force in the z direction is applied. The particles 
in the chain may circle the auxiliary particles, such that 
the couplings of the auxiliary particles to the particles 
in the chain equal. The proposal could also be realized 
with ultracold atoms in optical lattice by the method 
represented in [i^. We would like to notice that pre- 
venting any energy exchange between the chain and the 
auxiliary particles might be a challenging task in exper- 
iments. The requirement of uniform coupling also is ex- 
perimentally very challenging. Nevertheless, this paper 
established a link between the entanglement and critical- 
ity from the other aspect, which might shed light on the 
understanding of criticality from the auxiliary particles. 



and it may open up a new way to characterize and ex- 
perimentally detect quantum phase transition from the 
viewpoint of quantum information theory [2^, where two 
ancillas coupled to two particles at different sites in the 
chain are considered. By this model, the relation be- 
tween the state transfer quality and the spectral gap of 
the chain was established. 

To conclude, we have proposed a new method to study 
critical phenomena in many-body systems. The critical- 
ity was found to have a reflection in the entanglement 
of auxiliary particles that couple to it. Indeed, we have 
found that the entanglement change dramatically along 
the line of critical points. This dramatic change has been 
explained in terms of reflection of quantum phase tran- 
sitions, which lead to the coUectiveness in the couplings 
of the spin chain to the auxiliary particles. Quantum 
critical many-body system can also make a reflection for 
its criticality in decoherence of a quantum system which 
couples to itfiOl- In addition to shed light on the role of 
environment in quantum information processing, these 
pave a new way to study critical phenomena in many- 
body systems. The generalization of these results to 
a wide variety of critical phenomena and their relation 
to the critical exponents is a promising and challenging 
question which deserves extensive future investigation. 
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